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TIME-INHOMOGENEOUS JUMP PROCESSES AND VARIABLE 

ORDER OPERATORS 

ENZO ORSINGHER, COSTANTINO RICCIUTI, AND BRUNO TOALDO 


Abstract. In this paper we introduce non-decreasing jump processes with 
independent and time non-homogeneous increments. Although they are not 
Levy processes, they somehow generalize subordinators in the sense that their 
Laplace exponents are possibly different Bernstein functions for each time t. 
By means of these processes, a generalization of subordinate semigroups in the 
sense of Bochner is proposed. Because of time-inhomogeneity, two-parameter 
semigroups (propagators) arise and we provide a Phillips formula which leads 
to time dependent generators. The inverse processes are also investigated and 
the corresponding governing equations obtained in the form of generalized vari¬ 
able order fractional equations. An application to a generalized subordinate 
Brownian motion is also examined. 
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1. Introduction 

This paper is devoted to the study of non-negative, non-decreasing processes, 
say cr'^(f), with independent and non-stationary increments. We investigate their 
basic path and distributional properties with particular attention to the governing 
equations. Such processes are not Levy processes since the increments are not 
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stationary and therefore they consist in a generalization of subordinators. The 
relationship with Levy processes is highlighted by the fact that 


where A i— >■ n(A, t) is a Bernstein function for all fixed t > 0, i.e. 

nOO 

n{X,t) =Xb{t)+ / {l-e-^’^)(l>{ds,t) 

Jo 


( 1 . 1 ) 


( 1 . 2 ) 


under suitable assumptions on the couple {b{t),(j)(ds,t)). It turns out that such 
processes are stochastically continuous and a.s. right-continuous and therefore they 
are additive processes in the sense of [H, Definition 1.6]. We focus our attention 
on the case where 

i-t 

)ds 
iQ 


U{X,t) = [ /(A,s), 

Jo 

= J ^A&'(w)-|- J (l — z/(ds, w)^ dw (1.3) 

where A i—>■ /(A, s) is a Bernstein function for each fixed s > 0. Without the condi¬ 
tion of stationarity, the increments (T'^(t) —cr'^(s) have a distribution ^s,t depending 
on both s and t, and thus we have to deal with a two-parameter semigroup (also 
known as propagator). We propose a generalization of the Phillips’ formula for sub¬ 
ordinate semigroups in the sense of Bochner: if A generates a Co-semigroup of oper¬ 
ators Tt, t > 0, on some Hilbert space lo, then the map 1 1 —>■ Ts,tu = fg T^u ^B,tidw) 
is shown to solve the abstract Cauchy problem 


= -f 


(1.4) 


q{s) = u € Dom(H), 

where, for u S Dom(H), 

pOO 

—f{—A,t)u = b'{t)Au+ / {TgU — u)v{ds,t)- (1-5) 

Jo 

If /(x, t) = for a suitable function a{t) strictly between zero and one, we have 

-f{-A,t) = -{-Ar^*^ (1.6) 

and in this sense we say that the operator (11.51) is a variable order operator. Fur¬ 
thermore, under suitable assumptions, the one-dimensional marginal of cr'^(t) has 
a density /i(a;,t) satisfying equations of the form 


d d r _ 

—u{x,t) = J u{sA)v(.x - s,t)ds 


(1.7) 


where u(s,t) = z/((s, oo), t). In the case n(A,t) = fgX°‘^^'>ds, equation (11.711 be¬ 


comes 


d 


( 1 , 8 ) 


where the fractional derivative on the right-hand side must be understood in the 
Riemann-Liouville sense. The inverse processes 

L^{t) = inf {s > 0 : cr^(s) > t} 


(1.9) 
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are also investigated and the one-dimensional marginal of L^{t) is shown to solve 


dt 


— u{x,s)v{t - s,x)ds - Bt^xu{x, 


t) = --nix. 


t) 


( 1 . 10 ) 


where Bt^x is an integro-differential operator. Moreover, the composition of a homo¬ 
geneous Markov process via L^{t) is investigated and the corresponding governing 
equation is derived. 

There is a wide literature inspiring such results. First of all there are many papers 
devoted to the well-known natural relationship between stable subordinators and 
fractional equations. It is well known that the one-dimensional marginal of a stable 
subordinator cr“(t) solves the fractional equation dtu = —d^u. The density of 
the inverse of cr“ solves instead the time-fractional equation d^u = —dxU. The 
reader can consult 25, Chapter 8], [s^ and the references therein for further 
information on these topics. Note that equation (11.81) straigthforwardly generalizes 
these facts, and furthermore it turns out to be the governing equation of the so- 
called multistable subordinator studied by Molchanov and Ralchenko [s^. When 
the subordinator has a different Laplace exponent / then the governing equation 
can be written (see [^) as 

d d d r _ 

—u{x,t) = -b—u{x,t)-—J u{s,t)v{x - s)ds (l-H) 

and (HJD is the variable order generalization of ()1.11|1 since the kernel of the con¬ 
volution appearing in (O depends on the time variable. Concerning the inverse 
processes we have a classical result due to Baeumer and Meerschaert Q which states 
that the stochastic solution to a time-fractional Cauchy problem is a Levy motion 
time-changed via the inverse of an a-stable subordinator. This framework has been 
generalized to other subordinators and inverses by different authors. In Meerschaert 
and Schefher the authors pointed out the the scale limit of a continuous time 
random walk in is a Levy motion time-changed via the inverse of a general subor¬ 
dinator with Laplace exponent /(A). In this case is proved that the one dimensional 
marginal solves/(9t)u = Au where f {dt)u := C~^ [/(A)u(A) — A“^/(A)m( 0)] (t). In 
(i^ l the author studied operators of the form 




u{x, s)k{t 


s)ds 


( 1 . 12 ) 


under suitable assumptions on the Laplace transform C[k]. He relates (11.121) with 
the Poisson process time-changed via the inverse of a subordinator, described in 


In ^ it is shown that by means of the tail zz(s) of a Levy measure ;/(•) it is 


possible to write the equation of the inverse of a general subordinator as 

, d , , d 

b—u[x,t) J u{x,s)v{t-s)ds = -—u{x,t)- 


(1.13) 


The case in which D{s) is the tail of the Levy measure corresponding to the Laplace 
exponent /(A) = A^“ + 2/rA“, for a € (0,1/2) and /x > 0, is investigated in detail 
in [l^ and is related to finite velocity random motions. Note that equation (11.101) 
generalizes (11.131) in the sense the convolution-type operator on the left-hand side 
now depends on the space variable. We suggest for a nice summary of the 
classical theory of time-changed processes. For a different approach to time-changes 
and the corresponding fractional equations the reader can consult 
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An interesting particular case of the processes studied in the present paper is 
given by the multistable subordinator. Multistable processes provide models to 
study phenomena which locally look like stable Levy motions, but where the sta¬ 
bility index evolves in time. There are two different types of multistable processes 
(for a complete discussion see 27|). The first one is the so-called field-based process 
(see 0 ( 2 ^ 1 . which is neither a markovian nor a pure-jump process. The second 
one is the multistable process with independent increments (see 0 ) with Laplace 
exponent 


n(A,<) = [\^ 

Jo 




(1.14) 


which can be considered as the prototype of non-homogeneous subordinators. 
2. Non-homogeneous subordinators 


Our research concerns the cadlag processes 


a^{t) = b{t) -f 


E ' 

0<s<i 


is), 


t > 0 , 


( 2 . 1 ) 


where [0, 00 ) 9 1 1 —>■ b{t) is a non-negative, differentiable function such that 6(0) = 0, 
and e{s) is a Poisson point process in R+ with characteristic measure u{dx, dt). We 
will work throughout the whole paper under the following assumptions 

Al) v{ds,’) is absolutely continuous with respect to the Lebegue measure, i.e. 
there exists a density such that v{ds,dt) = v{ds,t)dt. Furthermore the 
family of measures {z/((is, is such that the function t 1 —>■ i'{ds,t) is 

continuous for each t. 

A2) for all t > 0, 

/ {x A 1 ) 1 '{dx, s)ds < 00 . (2.2) 

J (0,C!o) X [0,i] 

We call cr'^(t), t > 0, a non-homogeneous subordinator. 

Definition (12.111 consists in a slight generalization of the Levy-Ito decomposi¬ 
tion 0 which holds for non-decreasing Levy processes (subordinators). Therefore 
retains some important properties of the usual subordinators (that is the in¬ 
crements are independent and the sample paths are non-decreasing) but presents a 
fundamental difference consisting in the non-stationarity of the increments (whose 
distribution is here assumed to be time-dependent). Hence, the number of points 
of the poissonian process in any Borel set B C K."*" x K.+ of the form B = B x [s,t], 
where B C (0, 00 ), possesses a Poisson distribution with parameter 

m{B) = / v{dx,s)ds = / / i'{dx,w)dw. (2.3) 

In particular, the expected number of jumps of size [x, x dx) occurring up to an 
arbitrary instant t is given by 

4>{dx,t) = f iy{dx,T)dT. (2.4) 

Jo 

In view of (12.21) . which implies that 

pOO 

/ {x A l)(j){dx^t) <00 Vi > 0, (2-5) 

^0 
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we can apply Campbell theorem (see, for example, [23|, p. 28]) to the process (12.111 
in order to write that 


where 


Thus the function 


^OO 

n(A,f) =Xb{t)+ / {l-e-^^)(j){dx,t). 

Jo 

POO 

Ai-4n(A,t) = Xb{t) + / (^1 — cj){dx,t) 

Jo 


( 2 . 6 ) 

(2.7) 

( 2 . 8 ) 


is a Bernstein function for each value of f > 0. We recall that a Bernstein function 
/ is defined to be of class C°^ with (—l)"“^/("^(a:) > 0, for all n € N [^. Definition 

3.1] . Furthermore, a function / is a Bernstein function if and only if [39l . Theorem 

3.2] 

nOO 

f(X)=a + bX+ (1 - e"^") i/(ds) (2.9) 

Jo 

where a, 6 > 0 and iy{ds) is a measure on (0, cx)) such that 

poo 

/ {s A l)v{ds) < oo. (2-10) 

Jo 

Note that under Al) and A2), and the further assumption 

poo 

/ {x A l)i'{dx,t) < oo, Vt > 0, (2-11) 

Jo 

there exists a Bernstein function A i—?► f{X,t) such that (12.81) can be written as 

n(A,t) = J ^Xb\w)+ J [l — e~^^) i'{ds,w)^ dw = J f{X,w)dw. (2.12) 

In what follows the function s i—>• v{s,t) will denote the tail of the measure 
v{ds, t), i.e. 


7 ’^{'1, t) = oo), t), 7 > 0. 


(2.13) 


2.1. Paths properties. The process cr'^(t), t > 0, is the sum over a Poisson pro¬ 
cess, hence it has independent increments. As shown in the following theorems, 
(T^(t) is continuous a.s. and, under suitable conditions, strictly increasing on any 
finite interval. 

Theorem 2.1. The process is a.s. continuous, i.e., cr'^(t) = a.s., 

for each fixed t > 0. 

Proof. Observe that for h > 0 

Pr{lan(t)-an(t-h)l >e} =PrJ ^ g(s) > e 

I t—h<s<t 


<Pr< X! «(s)l{d-)<i} > I 


t —/t<S<t 
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+ Pri e(s)l{.(s)>i} > 


(2.14) 




Now since 


Pr< Y > I ^ 

—/i<s<t J t — h<is<t 

2 /■ /■' M 

= — / xv(dx,w)dw —;► U. 

^ JJ 0 

For the second term of (|2.14l) we have that 

Pr I ^ f I < 1 - e-^“ U-.., 


(2.15) 




'^0. 


(2.16) 


Continuity in probability implies that for any sequence f ^ it is true that there 
exists a subsequence such that cr^(tn) —>■ cr'^(t) a.s.. But since the processes cr'^(t) 
are cadlag the left limit must exist a.s. and therefore it must be equal to cr'^(t). 
Thus the theorem is proved. □ 


It is well-known that a Levy process is strictly increasing on any finite interval 
if the Levy measure is supported on (0, oo) (and hence is a subordinator) and has 
infinite mass, ;/(0, oo) = oo (see, for example, Theorem 21.3]). In our case a 
similar result is true. 


Proposition 2.2. Let W be a finite interval of [0,oo). If b'(t) > 0 for t £ W C 
[0,oo) the process (7^{t) is a.s. strictly increasing in W. Ifb'{t) = 0 for t £W but 
^((0,oo),f) = oo for all t £ W, then the process <J^{t) is a.s. strictly increasing on 

W. 


Proof. If b'{w) > 0 for w G [s,t] then it is clear that the process cr^(t) is strictly 
increasing. Now let b'{t) = 0. Observe that the Poisson point process {e{s),s > 0) 
is defined as the only point e(s) £ (0, oo) such that for a Poisson random measure 
(fi(’) on (0, oo) X [0,oo) with intensity v{dx,dt) it is true that 

¥^l(o,oo)x{t}(^*) “ ^{e{t),t)i.d,x) (2.17) 

where x £ (0,oo) x [0, oo). Fix an interval of time [s,f], then the probability that 
the process cr^(t), t > 0, does not increase in [s,t] is the probability that there are 
no points in (0,oo) x [s,t]. Let Ej, j G N, be a partition of (0,oo) with iy{Ej) < oo 
for all j £ N, then 


7^((0,oo) X [s,t]) = Y'i^i {{(0,00) nEj) x [s,t]) 


(2.18) 


and by the countable additivity Theorem [23|, p. 5] the sum (I2.18P diverges with 
probability one since = i^(0,oo) = oo. Since this is true for any interval 

[s,t] and if i^(0,oo) = oo the theorem is proved. □ 
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Remark 2.3. The process cr^{t), t > 0, defined as in (12.11) is an additive process 
in the sense of [s^, Definition 1.6]. This is because 

(1) it has independent increments (it is the sum over a Poisson process); 

(2) Xo = 0 a.s.; 

(3) it is stochastically continuous; 

(4) it is a.s. right-continuous with left-limit (Theorem [2T]). 

Stochastic continuity from the left (Item (3)) can be immediately verified as the 
stochastic continuity from the right (as in the proof of Theorem 12.II) . 


2.2. Time-inhomogeneous random sums. Since the compound Poisson process 
is the fundamental ingredient for the costruction of a standard subordinator, it is 
easy to imagine that random sums with time-dependent jumps play the same role 
in the definition of non-homogeneous subordinators. 

Let N{t), t > 0, be a non-homogeneous Poisson process with intensity g{t), t>0, 
and let Tj = inf{t > 0 : N{t) = j}. We consider the random sum 

Nit) 

Z{t) = ^ X{T,) (2.19) 

where X{Tj) is the positive-valued jump occurring at time Tj, having the condi¬ 
tional absolutely continuous distribution 

FT{X{Tj) G dx\ Tj =t} = tPidx,t), x>0, (2.20) 

with 

'tjj{dx,t) = I, Vt > 0. (2-21) 

The random variable Z(t) takes the value z = 0 with positive probability, and has 
a density for z > 0. Indeed 

Pr{Z(t) = 0} = Pr{7V(t) = 0} = e“ (2.22) 

and, for each z > 0 we have that 



Pr{Z(t) G dz} 


^ nt nt nt 

= ^/ / ... Pi{Z{t) e dz,Ti G dti,T2 e dt2,.. .Tn e dtn,N{t) = n} 

oo r^(t) 'I 

= /•••/ Pr< ^2 = ^2,... T„=t„,Af(t)=ni 

Jti Jt„^i [ 1=1 J 


X Pr {Ti G dti,T 2 G dt 2 , ... r„ G dtn, N{t) = n} 


°° j't j-t j-t ( " I 

/ / •■• / PJ'S ^dz} PrjTi G dti,...,Tn G dtn,N{t) = n} 

n=l -^0 dti J 

X{tj) G dz > Pr {Ti G dti,.. .Tn G dtn, N{t) = n} . (2.23) 


1=1 
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The first factor is given by the convolution integral 

Pr < > ^(tj) & dz > = dz / tj;{dxi,ti) ■ ■ ■ tp{dxm tn)d z — > Xj 
I I JiO.oo)" \ 




i=i 


(2.24) 


while the second one can be computed as follows 

Pr{Ti € dti,...Tn e dtn,N{t) = n} = g{ti)... g{tn)e~-^o ... dt^. 

(2.25) 

Then we have 

Pr{Z(f) G dz} 

= dz ^ / dti...dtn [ 'tjj{dxi,ti) ■ ■ ■1p{dXn,tn)s( Z - 

J[o,t]- JiO.oo)^ \ J 

xe-/o®M‘'"g(fi)...g(f„). 

We define the function 

= / g{T)i>ix,T)dT. 

Jo 

The density of Z{t) can be written as 

PT{Z{t) G dz} = dz e~-1^0 — W Xj \4){dxi,t) . . . <}){dxn,t) 

n=l -^(O.oo)" V J 

(2.28) 

and therefore its Laplace transform is given by 

]Eg-AZ(t) ^ pj.{z{t) =0}+ [ Pr{Z{t) G dz} 

J CO.oo') 


(2.26) 

(2.27) 


(2.29) 


where 


' (0,oo) 


e-^^Pr{Z(t) € dz} 


t ^ 1 f 

= g-Jo 9i-^)dT4){dxi,t) . . .4>{,dXn,t) 

J(o.oo)’* 


n—1 

CO 


= g-!o9(r)dr'^ }_ 


V-^(O.oo) J 


(2.30) 


= g- So 9i'^)d'^ e ^'‘<l>(dx,t) _ 

Combining all pieces together we have that 

gg-AZ(t) _ g-/o sMrf'r+/(o,oo) e-^^<p(dx,t) 

— g-So 9 (r)dr V’(da:u)+/(o.oo) So g(x)i/>(dx,T)dT 

= g-/(0,oo)(l-e“^*)<^'('^“.*) (2.31) 

with (j){x,t) of the form (j2.27l) . Observe therefore that in this case we have that 

i'(dx,t) = ip{dx,t)g(t) and u((0, oo), t) = g{t) < oo. (2.32) 
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2.3. Distributional properties. We remark that for each t > 0 the distribution 
of cr^(t) is infinitely divisible. In fact for each n G N, is given by 

the n-convolution of the probability measure of the r.v.’s associated to the Levy 
exponent 




nOO 

Jo 


^{ds, t) 


(2.33) 


where t is fixed. However, unlike what happens for Levy processes, such a distri¬ 
bution is not given by the n-th convolution of because the increments are 

not stationary. As pointed out in Remark [231 the process cr^(t) is an additive pro¬ 
cess. Therefore the fact that /xt(*) is infinitely divisible is also a consequence of [38, 
Theorem 9.1]. 

It is crucial to observe that cr^(t) can be approximated by means of a random 
sum of the form (I2.19L as stated in the following theorem. 

Theorem 2.4. Let be a non-homogeneous subordinator having Laplace ex¬ 

ponent 


poo 

U{X,t) = A 6 (t)+ / (1 

Jo 


— e 


— Xx 


)(j){dx,t) 


(2.34) 


as in () 2 . 8 I) and assume that s i—>■ v{s^t) is absolutely continuous on ( 0 ,oo) for all 
t > 0. Then there exists a process Z-^{t) of type \2.19\) such that, for j 0, we 
have 


b{f) Zj{t) —>• cr^(t). 

Proof. Let P be the function in (12.131) . Then 

, , vidx, t) , , 


(2.35) 


(2.36) 


is a probability distribution, because it is positive and integrates to 1. Let us 
consider the process 


Z^{t) = Y, X{T,) 
1=1 


(2.37) 


where Nj{t) is a non-homogeneous Poisson process with rate g^{t). We assume 
g-f{f) = b( 7 ,<) and use (12.361) to write 


Pr{X(Tj) € dx\Tj = t} = ipj{dx,t). 
In view of the discussion of Section 2.2 we have that 


(2.38) 


pz.,{dx,t) = (2.39) 


and 

So '>P-,idx,T)g-,{T)dT 

which converges to as 7 0 . 


(2.40) 

□ 
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Theorem 12.41 provides a method to construct an approximating process. As 
an example, let’s apply such a method to the multistable subordinator defined in 
Molchanov and Ralchenko 3^. In this case, the time-dependent Levy measure is 


(j){dx, t) = dx / 
Jo 


, , , ,,-dT (2.41) 

/o r(l - a(r)) 

for a suitable stability index r i—> Q;('r) with values in (0,1) in such a way that the 
conditions Al) and A2) are fulfilled. 

The Laplace transform thus reads 

= e“-^o (2.42) 

Being 

u( 7 ,t)=/ iy{dx,t) = — - — (2.43) 

J-r r(l-a(t)) 

the approximating process Zj(t) is based on a non-homogeneous Poisson process 


v-“(‘) 


with intensity g-y[t) = and jump distribution 


(2.44) 


A convenient way to deal with the non-homogeneity of the multistable process is 
to consider its localizability. We remind that (t) il localizable at t if the following 
limit holds in distribution (see, for example, 0 ))- 

+ (2.45) 


lim 

r—>-0 


rh(t) 


where Zt{T), T > 0 is the so-called local process (or tangent process) at time t. A 
fundamental property of Zt{T) is /i(t)-self-similarity, i.e. Zt{rT) = Zt{T) for 
r > 0. In the case where cr'^(t) is a multistable process, the local approximation 
at a fixed t > 0 is a stable subordinator with index a{t). By taking the Laplace 
tranform. 


lim Ee 
-- >-0 


= lim exp ■ 
1 —>-0 


= lim exp ■ 

r—>-0 


f‘t-\-rT 


rrA“(*) 


X 


^h{t) 


l(s) 


ds 


j{r) 


,-TA' 


,(t) 


(2.46) 

where the limit produces a non-trivial result by assuming that the similarity index 
is h{t) = l/a{t). 

Another way to approximate cr'^(t) is now given by means of stable processes. 
We split the interval [0, t] into n sub-intervals of length ^ and assume ai = a {^i)- 
We can write 






= lim e 

n—¥oo 




= lim I I ( 

n. —J" J- 




(2.47) 


and this proves that the following equality holds in distribution 

' i — 1 


^{t) = lim 

n—Voo • ^ 


2=1 


t] - (7^ 


-t 


= lim ( 

n—^oo < ^ 


2=1 


(2.48) 
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where cr“% 1 < i < n, are independent stable subordinators with index Ui = a (^*)- 
To conclude this section on distributional properties, we provide sufficient con¬ 
ditions for the absolute continuity with respect to the Lebesgue measure of the 
distribution of a non-homogeneous subordinator. 

Theorem 2.5. Let cr'^(t), t > 0 be a non-homogeneous subordinator with Laplace 
exponent 

pOC 

U{X,t) = Xb{t) + (l-e-^")0(ds,t) (2.49) 

Jo 

Suppose that s i—>■ v(s,t) is absolutely continuous on (0,oo) and furthermore assume 
Jp i^((0, oo), T)dr = oo for each t. The distribution of ^^{t) is absolutely continuous 
with respect to the Lebesgue measure. 

Proof. Consider the approximating process Z-fit) with distribution 

OO ^ 

pz,{dx,t) = + (2.50) 

n—1 

which converges weakly to the law of a^(t) as 7 —>■ 0 since Zjlt) o''^(t) as 
7 —>■ 0. Consider the Lebesgue decomposition of pz..,, written as 

pz-,{dx,t) = p%^{dx,t)-\-p%^{dx,t)-\-p^z.,idx,t). (2.51) 

By hypothesis we know that 

P%^{dx,t) = e-Jo^^-/’^)d^So{dx) (2.52) 

since s 1 —>■ z/(s, t) is absolutely continuous and therefore the Levy measure is abso¬ 
lutely continous with respect to the Lebesgue measure. Therefore by letting 7 —>• 0 
we observe 

{p% (dx, t) + (dx, t)) = e- Jo Hl,r)dr (2.53) 

which goes to zero as 7 —^ 0 since fgi7(0,T)dT = 00 and the continuity of the 
function 7 1 —^ f* 1 /( 7 , T)dT follows from the continuity of 7 1 —^ 1 ^( 7 , t) . □ 

2.4. The governing equations. Under the assumptions of the above theorem, we 
now derive the equation governing the density of a non-homogeneous subordinator. 
In the following, we will denote as q{x,t) the density of cr^(t), when it exists, i.e. 

Pt {a^{t) G dx} = q{x,t)dx. (2-54) 

Theorem 2.6. Let a^{t), t > 0, be a non-homogeneous subordinator and let the 
assumptions of Theorem 1 2. .51 hold. Then a Lebesgue density of (t) exists and 
solves the variable-order equation 

d d d r _ 

—g(x,t) = - 6 '(t)—g(a;,t) - — / q{s,t)v{x - s,t)ds, x>b(t),t>0, 

(2.55) 

provided that q{x, t) is differentiable with respect to x, subject to q{x, 0)dx = So{dx) 
for X >0, and q{b{t), t) = 0, for t > 0. 
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Proof. Now we consider the Laplace transform of the right-hand side of equation 
(j2.55|) and we get that 


d dr 

= Xb'{t)q{X,t) - b'{t)q{b(t),t) + XC[q*v\ (A) 


(A) 


= Xb'{t)q{X, t) + Xq{X, t) (A V(A, t) - b'it)) - b'{t)q{b{t),t) 


(2.56) 


where we used the fact that 

POO I 

J e~^‘'D{s,t)ds = jf{x,t) - b'{t). 
Therefore the solution to p.55l) has Laplace transform 
q{X,t) = 

which becomes, since q{b{t),t) = 0, 
q{X,t) = 

and coincides with 


(2.57) 


(2.58) 

(2.59) 
□ 


If cr^{t) is a multistable subordinator with index Qf(<), we have 


v{x,t) = 


r-a(t) 


and the governing equation reads 


^ r(l — ait)) dx Jq 


r(l-a(t))’ 

d 




1 


:dy. 


(2.60) 


(2.61) 


{x — y)“d) 

Keeping in mind the definition of the Riemann-Liouville fractional derivative of 
order a G (0,1) 


5“ 1 d r Mfj), 

r{l-a)dxJo {x-y)^‘^^' 


we can write ()2.61[l as 
d 


ga(t) 


0 < aft) < 1, a: > 0, 


(2.62) 


(2.63) 


where g^dt) is the Riemann-Liouville derivative of time-varying order aft). Then, 
by taking inspiration from 4^ Definition 2.1], we define the generalized Riemann- 
Liouville derivative with kernel u(x, t) as 


V^ft) q{x,t) = -^ q{s,t)v{x - s,t)ds (2.64) 

where the operator ^^{t) acts on the variable x but also depends on t. Using this 
notation, we say that the density of a non-homogeneous subordinator solves the 
following Cauchy problem 


■§tq{x, t) = - V^{t) q{x, t), t>0, 

q(x,0) = S(x). 


(2.65) 
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It is useful to define also a generalization of the Caputo fractional derivative as 

(t) = J ■^qis,t)v(x - s,t)ds. (2.66) 

If X !->■ q{x, t) is absolutely continuous on [0, oo) then V'^ (t) exists a.e. for all t > 0, 
and the following relationship holds 

t) = 9(0, t)V{x, t) + {t) q{x, t) (2.67) 

whose proof can follow [dll Proposition 2.7]. Formula (12.671) is a generalization 
of the well-known classical relationship between Caputo and Riemann-Liouville 
derivatives [l^, page 91]. 


3. The inverse process 


In this section we consider the process 

L^{t) = inf {x > 0 : <J^{x) > t} (3.1) 

where <t^{x) is a non-homogeneous subordinator without drift, namely b'{x) = 0 
for all X. We throughout assume that 

p((0, oo),t) = oo for alH > 0 (3.2) 

and that 


s I—>■ iy{s,t) = p((s, oo),t) is an absolutely continuous function on(0, oo). (3.3) 


By using Theorem [2T] and Remark 1 2.2 1 it is clear that the process is well defined 
as the inverse process of cr^it). Observe that, a.s., L^{a^{t)) = t since L^{a^{t)) = 
inf {s > 0 : cr^(s) > cr^(t)} and, under (|3.2I1 . the process cr^(f) is strictly increasing 
on any finite time interval fProposition 12.21) . In the following, we denote by a: i—>■ 
l{x,t) the Lebesgue density of when such a density exists. The inverse of a 

classical subordinator has a Lebesgue density ([i^ Theorem 3.1]). We provide here 
an equivalent version of Theorem 3.1] valid for non-homogeneous subordinators. 

Theorem 3.1. Under the assumptions (lO) and (1331) the process t > 0, 

has a Lebesgue density which can be written as 

X 1-^ l{x,t) = / q{s,x)D{t — s,x)ds, forallt>0. (3-4) 

Jo 


Proof. Define 


and 


L{z,t) = / l{x,t)dx 
Jo 


R{z,t) = PT{L^{t) < z} . 


(3.5) 

(3.6) 


We will show that L{z,t) = R{z,t). By using the convolution theorem for the 
Laplace transform we have that 

£(^,A) = (3.7) 

The use of the relationship 

Pr {L^it) > a;} = Pr {cr"(x) < t} 


(3.8) 
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leads to the Laplace transform 

pOO pCO 

/ e~^*R{z,t)dt = / e~[l — Pr {cr^{x) < t}) dt 
Jo Jo 


1 1 

“ A “ A® 


-n(A,a:) 


(3.9) 

(3.10) 


Therefore we have proved that 

nOO nOC 

/ e~^*L{z,t)dt = / e~^*R{z,t)dt 
Jo Jo 

and thus in any point of continuity it is true that 

i?(z,t) = L{z,t). (3.11) 

If we prove that 1 1 —>■ R{z, t) and 1 1 —>■ L{z, t) are continuous functions then we have 
proved the Theorem for all t. Note that under (13.21) the process <T^{t) is strictly 
increasing on any finite interval in view of Proposition 12.21 Therefore the process 
L^{t) is a.s. continuous and therefore it is also continuous in distribution. This 
implies that t i-A R{z, t) is continuous. Now we show that 1 1 —>■ L{z, t) is continuous. 
Note that, for h > 0, 

\l{x, t + h) — l{x,t)\ 

pt-\-hj pt 

/ q{s^x)9{t-\-h — s^x)ds — / q{s^x)9{t — s^x)ds 

Jo Jo 

pt pt-\-h 

/ q{s,x) {D{t + h — s, x) — D{t — s,x)) ds + / q{s,x)D{t + h — s,x)ds 

pt pt-\-h 

< / q{s,x) \D{t + h — s,x) — D{t — s,x)\ds + / q{s,x)p{t + h — s,x)ds 


pZ pt-\-h 

= / q{s^x) {D{t — s,x) — + h — s,x)) ds / q{s,x)D{t + h — s,x)ds. 

(3.12) 

Since under (1531) the function s i-A u(s, •) is absolutely continuous and since 

v{t — s,x) — 9(1 — s + h,x) < v(t — s,x) (3.13) 

and 

/ v{s,x)ds<oo, (3.14) 

Jo 

the first integral in (I3.12|) goes to zero by an application of the dominated conver¬ 
gence theorem. The second integral is for any oo > z > t and sufficiently small 
h 

/ t-\-h pz 

q{s,x)D{t + h — s,x)ds = / q{s,x)t(^tt+h){s)Pit + h — s,x)ds. (3.15) 
Now since 

+ h - s,x)ds < (?(s,x)l(t_ 2 )(s)b(t - s,x)ds (3.16) 

and 


q{s,x)l(^t,z){s)i^{t — s,x)ds < oo 


( 3 . 17 ) 
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another application of the dominated convergence theorem shows that the second 
integral in (j3.12|) goes to zero. For h < 0 the arguments are similar. This completes 
the proof. □ 

Theorem 3.2. If x ^ v{t,x) is differentiable and if the density x i— l{x,t) is 
differentiable then l{x,t) solves the equation 
d 

—l(x,t)=S(x)T'(t,x)—'D?(x)l(x,t) — Btxl(x,t), x>0, (3.18) 

ox 

in the sense of distributions, namely it solves pointwise the Cauchy problem 


-§^l{x, f) = -Vf-{x) l{x, t) - Bt,xl{x, t) a; > 0 


(3.19) 


\^l{0, t) = iy{t, 0) 

where 2?^(x) is the generalized Riemann-Liouville derivative acting on t and de¬ 
pending on X > 0, and Bt^x is an operator acting on both t and x defined as 

C d dr 

Bt,xl{x,t)= J ds—iz{t-s,x)— J l{x',s)dx'. (3.20) 

Proof. We can adapt [2^, Theorem 8.4.1] to our case. It is sufficient to derive both 
sides of (lOl) and apply 


d . d , . 


(3.21) 


to obtain 

A 

dx 


f d _ f d _ 

l{x,t) = / -—q{s,x)n{t—s,x)ds-\- / q{s,x)-—n(t — s,x)ds 
Jo iJx Jq dx 

r d _ r d _ dr 

= — -^l(x, s)iy{t — s,x)ds — / ds-^ult — S.x)-;:— / l{x',s)dx' 

Jo ds Jo dx ds Jq 

C Q d 

=—V?ix) l{x,f) — / ds-—v(t—s,x)-— / l(x',s)dx' 

Jo Jq 

r d dr 

= 5{x)v{t,x)—Vf'ix) llx,t) — / ds-^Vit — s,x)-— / l{x',s)dx' 

Jq dx ds Jq 


where in the last step we referred to (12.671) . 


(3.22) 

□ 


Remark 3.3. Non stationarity is here expressed by the term Bt^xl{x,t), which 
vanishes in the case of the inverse of a classical subordinator, and by the fact that 
the kernel of V^ix) depends on x > 0. 

In the case of the inverse of a classical stable subordinator, Theorem l3.2l obviouslv 
leads to the well-known Cauchy problem eq (5.7)] 


dx 


l{x,t) = —^l{x,f), X > 0,t > 0 


(3.23) 


— r(l-a)’ 

for a € (0,1). 

We remark that a first study on time fractional equations with state-dependent 
index appears in 0 where the following equation is considered 

da{k) 

= -0 (j)k{t) - Pk-l{t)) . 


(3.24) 
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The fractional derivative is meant in the Dzerbayshan-Caputo sense and the second 
member arises by writing the forward equations with the adjoint of the generator 
of the Poisson semigroup. 

3.1. Time changed Markov processes via the inverse of non-homogeneons 
subordinators. We now consider the composition of a Markov process with the 
inverse of a non-homogeneous subordinator. Let X(u), u > 0 he a Markov process 
in such that X(0) = y a.s. and 

Pr{X(u) € dx} = p(x,y,u)dx. (3.25) 

We assume that p(x, y, u) is a smooth probability density satisfying the following 
Cauchy problem: 

\^P = u>0, 

\p{x,y,0) = S{x - y), 


(3.26) 


where Sx is the adjoint of the Markovian generator acting on the variable x. More¬ 
over let t > 0, be the inverse of a non-homogeneons subordinator, with 

density as in Theorem 13.II 

Pr{T^(t) S dx} = l{x,t)dx. (3.27) 

By assuming that X{t) and L^{t) are independent, we study the composition 
X{L^{t)), having distribution 

poo 

Pr{XG dx} = / Pt{X{u) G dx}Pi{L^{t) G du}. (3.28) 

Jo 

Then X{L^{t)) has a smooth density, defined as 


p<X> 

gix,y,t)= p{x,y,u)l{u,t)du. 
Jo 


(3.29) 


By using simple arguments, we now derive the governing equation for ()3.29p . 

Proposition 3.4. Under the above assumptions, the density i3. 29\) . for t > 0, 
solves the following equation in the sense of distributions: 

poo 


(«) [p{x,y,u)l{u,t)\du =5{y - x)v(t,0) + Sxg{x,y,t) 


f 


p{x,y,u)Bt^ul{u,t)du. 


(3.30) 


Proof We have 

poo poo 

/ 'D^{u) \p{x,y,u)l(u,t)]du= / p{x,y,u)'Df'{u) l{u,t)du (3.31) 

Jo Jo 

and by using (13.191) and (I3.26|) . which hold for positive times, we can write 

pOO 

/ [p{x,y,u)l{u,t)\du 

Jo 

/.oo Q noo 

= — lim / p(X, y,u)—l{u,t)du — \im / p{x,y,u)Bt ul(u,t)du 
e^oj^ du ^^OJe 

d 

= — lim [p(a:, y, m)1(u, t)]“-I-lim / —plx,y,u)Uu,t)du 
e->o e-j-o au 
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pOO 

- lim / p{x,y,u)Bt^ul{u,t)du 

pOO pOO 

= 5{y - x)V{t,0) + \im Sxp{x,y,u) l{u,t)du - / p[x,y,u)BtJ{u,t)d' 

Jo 

poo 

= S{y - x)v{t, 0) + Sxg{x, y,t)- / p{x, y, u) Bt^u l{u, t)du 

Jo 

and the proof is complete. 


(3.32) 

□ 


Remark 3.5. In the case X{t) is a Brownian motion starting from y and L^{t) is 
the inverse of a multistable subordinator with index a{x) € (0,1) we have x) = 
r{i-a(x)) ■: governing equation reads 

OO Qa{u) 

0 ^^ [p{x, y, u) l{u, t)] du = -A^g{x, y, t) + S{y - 

f°° 1 (y-^d 

e Bt^ul{u,t)du x>0 


L 


\/2'ku 


(3.33) 


where a(0) = oq and 

[* \ d _ dr 

Bt^uKx,t) = / ds —v{t - s,u)— / l{w,s)dv 

Jo \_OU OS Jq 


= ds 


d (t-s)-“(“) d 


l{w, s)dw 


(3.34) 


Iq idu r{l — a{u)) ds Jq 

Note that (13.3311 is a generalization of the well-known fractional diffusion equation 
to which it reduces when u i—>■ a{u) is constant, that is 


d° 


r 


:9 


-5{x -y) = A^g 


(3.35) 


dfa^ r(l — a) 

Remark 3.6. Consider the case where the Markov process is a deterministic time, 
namely the starting point is y = 0 and X{t) = t. In this case we have Sx = so 
that the governing equation becomes, for a; > 0, 

/ 'Dj^{u) [p{x,y,u) l(u,t)]du = S{x)V{t,0) — —g — / (5(a; — Z(it, t)(itt 

Jo dx Jq 

(3.36) 

and obviously coincides with that of since the probability density of X{u) is 
p{x, 0, u) = S(x — u). 


4. Non-homogeneous Bochner subordination 

We consider in this section a generalization of the Bochner subordination. We 
recall here some basic facts. Let Tt be a Cp-semigroup of operators (the reader can 
consult for classical information on this topic) i.e. a family of linear operators 
on a Banach space (f8, IHI^) such that, for all m G IB, 

(1) Tq-u = u 

(2) TtT^u = Tt+sU, s,t > 0, 

(3) limt_>o llTtit - = 0. 
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Let (^,Dom(zl)) be the generator of Tt, i.e. the operator 


Au := lim 

t —^0 


Ttu — u 
t 


defined on 


Dom(A) = 


M G $ 


lim 

t —^0 


TtU — u 
t 


exists as strong limit 


(4.1) 


(4.2) 


and let ||rtu||^ < 

Let be a convolution semigroup of sub-probability measures associated with 
a subordinator, i.e. a family of measures {fJ-t}t>o satisfying 

(1) //((O, oo) < 1, for alH > 0, 

(2) fit * fig = /it + S, 

(3) limt^o = <5o vaguely, 
and such that 


£N(A) = 


(4.3) 


where 


poo 

/(A) = a + bX + (l — v{ds) 

Jo 

is a Bernstein function. The operator defined by the Bochner integral 


r/u = 


TgUfitids), 


M G 


(4.4) 


(4.5) 


is said to be a subordinate semigroup in the sense of Bochner. A classical result 
due to Phillips 371 states that t/ is again a Co-semigroup and is generated by 


—f{—A)u:=—au + bAu+ / {TgU — u) i'{ds) 


(4.6) 


which is always defined at least on Dom(A) Theorem 12.6]. 

In order to extend such a result to non-homogeneous evolutions, a generalization 
of the notion of one-parameter semigroup is needed. Let (IB, ||.||(g) be a Banach 
space. A family of mappings Ts,t from 18 to itself, defined by the pair of numbers s 
and t (such that 0 < s < t), is said to be a propagator (two-parameter semigroup) 
if for each m G 18, 25, Section 1.9] 


(1) Tt,tu = M, for each t > 0; 

(2) Ts,tTr,sU = Tr,tU, iOT r < S < t] 

(3) lim.s^o\\Ts+s,tu-Ts,tu\\^ = lim5_),o ||7^,t-i-5M - 7^.tu||<8 = 0; 

It is obvious that a propagator Ts,t reduces to a classical one-parameter semigroup 
in the case where it only depends on the difference t — s. 

Let cr^(t), t > 0, be a non-homogeneous subordinator and consider the measures 
fis,t{-) corresponding to the distribution of the increments cr'^(t) — o''^(s) which are 
obviously such that 






(4.7) 


as can be ascertained by applying the Campbell theorem to a^{t) — cr^(s) under 
the assumption (12.111) . Therefore, it is easy to verify that the family of measures 
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{/is.i(*)}o<s<t forms a two-parameter convolution semigroup of probability mea¬ 
sures since, from the independence of the increments and (14.7L we get * IJ>r,s = 
r < s <t. Consider the operator defined by the Bochner integral on 05 


nOO 

Ts,tU= / Hs,tidu)- 

Jo 


(4.8) 


The family of operators {7s,t}o<s<t forms a two-parameter semigroup of operators 
on 05, i.e., ()4.8(1 is a propagator. This can be easily ascertained by observing that 
for all M £ 05 






He,tidw) 


/ Tw'U ^jLr,s{dw') 

Jo 

Tw-\-w'^ l^r,si^dw ^ 

TpU )J.r,sid{p - w)) fj.s,tid'w) 
TpU f fis,tidip - w)) Hr,sidw) 

Jo 

) 

TpU 


OO poo 


0 ^0 
OO poo 


0 J w 

OO pp 


— Tr.tU. 


(4.9) 


We consider here the case where the generator {A, Dom(yl)) of T* is a self- 
adjoint, dissipative operator on an Hilbert space (f), (v^) and thus we have that 
117)^11^5 < llwll^j (see, for example, [l^, Section 2.7] and [H, Chapter 11] for classical 
information on linear operators on Hilbert spaces). Recall that an operator is said 
to be dissipative if {Au, u) < 0 for all u £ Dom(H) and 

Dom(H) = |it £ < ooj . (4-10) 


Theorem 4.1. Let the above assumptions (including (12.111) ) be fulfilled. The family 
of operators Ts.t acting on an element u € is a bounded propagator on Sj and for 
u £ Dom{A), the map 1 1 —>■ Ts,tu solves, 


0<s<t, 

|( 7 (s) = M £ Dom{A), 

where the family of generators {—/(—^, t)}(>Q, can be defined as 

pOO 

-f{-A,t)q ■= b'{t)Aq+ / [T^q- q)v{ds,t), 

Jo 

a Bochner integral on Dom{A). 


(4.11) 


(4.12) 


Proof. First note that 

poo 

< / \\Twu\\f^Ps,tidw) <\\u\\^, 
Jo 

and therefore is bounded. 


(4.13) 
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Then we recall ([s^, Theorem 11.4] and [l^ Theorem 2.7.30]) that within such 
a framework we have by the spectral Theorem that for u € Dom(A) 


Au = / XE{dX)u 

J ( —C!O,0] 


(4.14) 


where E{B) : Dom(A) Dom(A), B a Borel set of K., is an orthogonal projection¬ 
valued measure supported on the spectrum of A defined as 


E{B)u := j E{dX)u. 
Jb 


(4.15) 


Therefore since from (14.1411 it is true that for a function $ : (—oo, 0] i—?► R 

<^{A)u = f ‘^>{X)E{dX)u (4.16) 

J f —oo.Ol 


we have that 


Ttu = / e*^E{dX)u. 

J ( — 00 , 0 ] 


(4.17) 


We now verify that Ts,tTr,sU = 'Tr,tu, r < s < t, in order to show that for 
all u G the operator Es,t is a propagator since the other defining properties as 
trivially verified. We have that, for all u G ij, 


^oo noo 




0 ^0 
oo poo 


T^+pU {dp) 

E{dX)E{dQ)u pr,s {dw)ps,t {dp) 


0 Jo J { — oo,0] J { — oo,0] 


Sr f{->^,'w)dw - /* f{-Q,w)dt 


'E{dX)E{dg)u 


( — 00 , 0 ] J ( — 00 , 0 ] 


= -/r f{->^,w)du 


( — 00 , 0 ] 


e ■—'-"E{dX)u 


Hr,t{d'w)E{dX)u 


Jo J ( — 00,0] 

pOO 

= / Tu,U PrAd't^) 

Jo 

— 

For a function u such that 

f \f{-X,t)f{E{dX)u,u)<oo 
J ( — 00,0] 

the representation (14.121) can be shown to be true: use (14.161) to write 

-f{-A,t)u = -[ f{-X,t)E{dX)u 

J ( — 00 , 0 ] 

= — f (—b'{t)X+ f (l — v{ds,t)] E{dX)u 

= [ b'{t)XE{dX)u+ [ [ -1) E{dX)ui'{ds,t) 

J(-oo.O] Jo J(-oc.O] 


(4.18) 

(4.19) 
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(4.20) 


pOO 

= b'{t)Au+ / {TgU — u) ^{dsA)- 
Jo 

Now we show that (I4.20p is true for any u € Dom(A) 

pOC 

\\fi-At)u\\^ <b'{t)\\Au\\^ + / ||T,m-u||^ j/(ds,t) 

^0 

poo 

<b'{t)\\Au\\p^+J s\\Au\\p^iy{ds,t)+ 2 J ||u||^^ <). (4.21) 


Now note that 


poo 

Ts,tU = / T^u pLs^t{dw) 

Jo 

= [ f e'^^E{d\)u 

Jo J(-oo,0] 


/is.t(dw) 


/ 


■ !l f{-^,T)dT , 


( —C3O,0] 


g jsjv E{dX)u 


(4.22) 


where we used ( 4.161) . The fact that Ts,t maps Dom(A) into itself can be ascertained 
by using again 39|, Theorem 11.4] for saying that E{’) maps Dom(A) into itself and 
furthermore, since E{I)E[J) = E{I fl J) for any /, J Borel sets of K., we observe 
that for any u G Dom(A) 


Ts 


tAu = [ e-J‘f^-^''"^‘^^E{dX) [ nE{dn)i 

J ( — 00 , 0 ] J ( — 00 , 0 ] 

= f Xe-JJ^^-^’^^‘^'^E{dX)u 

J ( —oo.Ol 


J ( — 00 , 0 ] 

= ATs,tu. 


fiE{dfi) 




( — 00)0] 


e '"-'-^E{dX)u 


Now note that the equality 
d 


dt 


Ts,tu = -f{-A,t)Ts,tu, 0 < s < t, 


(4.23) 


(4.24) 


must be true in the sense of (14.161) and indeed by using (14.221) we have that, for 
u £ Dom(A), 


jLp- !l 


( — 00 , 0 ] 


dt 


[ f{-X,t)e-J‘^^-^’^^‘^'^E{dX)i 

J ( — 00 , 0 ] 


d _ 


f{-d,t)E{d^i) 

f ( — 00 , 0 ] 

= - f{-AA)rs,tu 

where we used again [s^. Theorem 11.4]. 


E{dX)u 


g-Js /(-A.u))dti 


( — 00 , 0 ] 


E{dX)i 


(4.25) 

□ 
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4.1. Time-changed Brownian motion via non-homogeneous subordina- 
tors. In this section we provide some basic facts concerning Brownian motion 
time-changed with a non-homogeneous subordinator. This is the immediate gen¬ 
eralization of the classical subordinate Brownian motion: the reader can consult 
[3 HH; [12 Hi for recent developments on this point. Therefore we assume now that 

Ttu = W.=^u{B{t)), t>0, (4.26) 

for u S where B is an n-dimensional Brownian motion starting from x S M”. 

We have therefore the formal representation 

Ttu = (4.27) 


where A is the n-dimensional Laplace operator such that 

Dom(A) = |n G L^(]R") : ||Au||j^ 2 (Rn) < oo| . (4.28) 

Therefore we get 

Ts^tu = E^u{B{a^{t)-a^{s))) , 0 < s < t,it G (K”). (4.29) 

Consider, for example, the case of a multistable subordinator, where n(A,t) = 
f* fQj. suitable choice of a(s) with values in (0,1). Then A >->■ A“(®^ is a 

Bernstein function for each s > 0, and Theorem 14.II leads to 

- (-A)“(*) u = r (TsU - u) (4.30) 

r(l -a(t)) Jo 

for a function u G Dom(A). Note that in this case we have a Brownian motion 
composed with the multistable subordinator whose increments have characteristic 
function 

]Egj?S(iTn(t)-(T“(s)) ^ ( 4 . 31 ) 

By following, for example, [HI Section 3.1] the generator ()4.3Q|1 can be also defined 
as 

- (-A)“(‘) u = ll^f (4.32) 

with 

Dom ^(—A)“^*^^ = G (R") : J < oo, for each t > 0 

(4.33) 


In general, we observe that for any non-homogeneous subordinator we can write 

]Eg»«s(<^“(i)-'^”(s)) = (4 34) 


and we can adapt [l2, Example 4.1.30] to write 


“/ ( I 11 = - 


1 






,t]u{C)d^ (4.35) 


with 




,tj u{^)d^ < oo,Vt > 0j . 

(4.36) 
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Therefore, we have by Theorem 14.11 the structure of the solution to a sort of 
diffusion equation 




q{t). 


(4.37) 


Remark 4.2. These last remarks have some relationships with Hoh’s s ymb olic 
calculus (discussed in [l^, 17, Chapter 6 and 7]) and in particular with (1^ . In 


full generality one can consider a Fourier symbol of the form 
where / : [0, oo) x R" i—>■ R. is a Bernstein function for each fixed x G 


(4.38) 
C and 


12 


t(x,^) is a symbol in the Hoh’s class. Then under some technical assumptions 
Theorem 2.4] it is true that the (variable order) pseudodifferential operator 

Q(")u(x) := (27r)-”/2 / e^-<f{t{x,0,x)d^ (4.39) 

generates a Feller semigroup (on Cq (R")). The Fourier symbol of our processes 
can be of the form /(t(x,^),t) but the dependence on t clearly originates non- 
homogeneous Markov processes and therefore two-parameter semigroups. 

We investigate here the mean square displacement i.e. the quantity 

^{t) = [ ||x||^ Pr {R (cr^(t)) € dx} . 

jR" 


(4.40) 


Roughly speaking, a stochastic process is said to have a diffusive asymptotic be¬ 
haviour when fM{t) ^ Ct i.e. the mean square displacement grows linearly with 
time. When M{t) ^ t°‘, a G (0,1), the process is said to be subdiffusive, while 
> 1 it is super-diffusive (the reader can consult 


if a 


for an overview on 


anomalous diffusive behaviours). Here it is interesting to note that the mean value 
of the Levy measure, namely wi/(dw, t) determines under which conditions the 
asymptotic behavior is respectively diffusive, sub-diffusive or super-diffusive. 


Proposition 4.3. We have the following behaviours. 
(1) If and only if 


/ OO 

w<f>{dw, t) 


< oo for 0 < t < to < oo 


it is true that M{t) < oo for all t < to 

(2) Under 1)4.411) for to = oo, we have that 

fW(t) . 

0 < lim - = C < oo if and only if lim / wv(dw, t) = C 

t—)-oo t t—^OO Jq 

Mit) . 

lim - = oo if and only if lim / wv[dw, t) = oo 

t—¥00 t t—¥00 Jq 

(3) Under (14.411) for to = oo, if 


then 


pOO 

lim / wv{dw, t) = 0 
Jo 


lim = 0. 

t—^OO t 


(4.41) 

(4.42) 

(4.43) 

(4.44) 

(4.45) 
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Proof. Observe that under (I4.41|l 

iW(t) = f ||x||^ / Pr {i?(s) G dx} Pr{(T'^(t) G ds} 

dR" Jo 

/'CO 

= n sPr S ds} 

Jo 


- - n— 


A=0 


POO 

= n w(j){dw,t) < oo if t < to- 
do 


(4.46) 


Observe that the last integral in (14.4611 converges only under (I4.41|) . Now note that 

M{t) n 'w(j){dw,t) 


lim 

i—)-oo t 


= lim 

t—^co 


n fg wv(dw, s)ds 


and therefore the proof of Item (2) and (3) is easy to be done. 


(4.47) 

□ 


A time-change by means of a multistable subordinator leads in this case to a 
process with 9v[{t) = oo for any t as a consequence of Item [T] of Proposition 14.31 
Consider now the measure 


v{ds,t) = s ds 


(4.48) 


for a function a{t) > 0 such that Al) and A2) are fulfilled. The associated Bernstein 
functions become, for each t > 0, 


/(A,t) = log 1 + 


A 

a{t) 

and in view (14.481) we can compute 

, . /"* dr 

=n —r^. 

Jo a{T) 

Observe that Proposition 1131 leads to the study of the limit 

poo 


f 

lim / 

t^oojg 




dw = lim 


I 


t->-oo a(t) 


(4.49) 


(4.50) 


(4.51) 


therefore the asymptotic behaviour of M (t) in this case depends on the asymptotic 
behaviour of aft). 

If instead, for functions a(t) strictly between zero and one and 0(t) > 0 as in 
AI) and A2), 


v{ds,t) = 




ds 


(4.52) 


r(l - aft)) 

then the Bernstein functions are a generalization of the Laplace exponent of the 
relativistic stable subordinator 

/(A, t) = {X + 6»(t))“^‘^ - 6»(t)“(‘) (4.53) 

and the asymptotic behaviour of the fMft) is determined in this case by the limit 


....X , -ds = lim ait) 9{t)°‘^*^^ 

t^oo jg r(l - aft)) t^oo ^ ^ 


lim 


(4.54) 
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The explicit form of M{t) is here 

M{t) = n [ f \ . dr. (4.55) 

Remark 4.4. For the reader who is familiar with fractional (anomalous) diffusion 
processes we give here some intuitive flashes of insight. The subordinate Brownian 
motion, as well as its generalization proposed above is, in general, a non pathwise 
continuous process. Hence it is not a diffusion. However the subordinate Brownian 
motion is sometimes included in the class of the so-called fractional diffusions (or 
anomalous diffusions) since the discontinuity of the sample paths is introduced via 
a time-change and therefore it develops in an operational time. If, now, we denote 
such an operational time as t* then the generalized subordinate Brownian motion is 
the process B ((T^(t*)) where the external time is t = Therefore it must be 

clear that the parameters must depend on the operational time t* and not on the 
externally measured time t. In the fractional case this dependece must be meant 
as t* !->■ a{t*) as well as in the tempered case we have t* ^ 
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